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We derive the effective equations for the out of equilibrium time evolution of the 
order parameter and the fluctuations of a scalar field theory in spatially flat FRW 
cosmologies. After setting the problem in general we propose a non-perturbative, 
self-consistent Hartree approximation. The method consists of evolving an initial 
functional thermal density matrix in time and is suitable for studying phase transi- 
tions out of equilibrium. The renormalization aspects are studied in detail and we 
find that the counterterms depend on the initial state. We investigate the high tem- 
perature expansion and show that it breaks down at long times. The infinite time 
limit is computed for de Sitter spacetime. We obtain the time evolution of the initial 
Boltzmann distribution functions, and argue that in the Hartree approximation, the 
time evolved state is a "squeezed" state. We illustrate the departure from thermal 
equilibrium by studying the case of a free field in de Sitter and radiation dominated 
cosmologies. It is found that a suitably defined non-equilibrium entropy per mode 
increases linearly with comoving time in a de Sitter cosmology, whereas it is not a 
monotonically increasing function in the radiation dominated case. This work has 
been done in collaboration with D. Boyanovsky and R. Holman. 

1 Introduction 

The usual assumption in inflation |I[ is that the dynamics of the spatial zero mode 
of the (so-called) inflaton field is governed by some approximation by the effective 



potential which incorporates the effects of quantum fluctuations of the field. Thus 
the equation of motion is usually an ordinary non-linear differential equation of the 
form: 

+ 3^ + Kff(0)=O (1) 

where a(t) is the expansion factor in the metric: 

ds 2 = dt 2 -a 2 (t) dx 2 (2) 

The problem here is that the effective potential is really only suited for analyzing static 
situations; it is the effective action evaluated for a field configuration that is constant 
in time [0. Thus, it is inconsistent to use the effective potential in a dynamical 
situation. Notice that such inconsistency appears for any inflationary scenario (old, 
new, chaotic, ...). 

More generally, the standard methods of high temperature field theory are based 
on an equilibrium formalism]^, |J; there is no time evolution in such a situation. Such 
techniques preclude us from treating non-equilibrium situations such as surely exist 
for very weakly coupled theories in the early universe. 

More precisely, the effective potential develops an imaginary part around points 
<f)Q where V"((po) < [<f>o : false vacuum]. This fact leads to long- wavelength modes 
with imaginary frequencies uj{k) = \Jk 2 + V"((po) for < k < \J\V"((po)\. When the 
field starts localized around 0o , these long-wavelength unstable modes push the 
system to evolve on time away from O - The smaller is the non- linear coupling, the 
better is this small fluctuations approximations and the more important are these 
unstabilities. 

In collaboration with Daniel Boyanovsky and Rich Holman we try to rectify this 
situation by addressing three issues^, ||, [/J ||: a) obtaining the evolution equations 
for the order parameter including the quantum fluctuations, b) studying departures 
from thermal equilibrium if the initial state is specified as a thermal ensemble, c) 
understanding the renormalization aspects and the validity of the high temperature 
expansion. 

Our ultimate goal is to study the dynamics of phase transitions in the early uni- 
verse, in particular, the formation and evolution of correlated domains and symmetry 
breaking in an expanding universe. From some of our previous studies on the dynam- 
ics of phase transitions || |(| in Minkowski space, we have learned that the familiar 
picture of "rolling" is drastically modified when the fluctuations are taken into ac- 
count. As the phase transition proceeds fluctuations become large and correlated 
regions (domains) begin to grow. This enhancement of the fluctuations modifies sub- 
stantially the evolution equation of the order parameter. Thus the time dependence 
of the order parameter is not enough to understand the dynamical aspects of the 
phase transition; it must be studied in conjunction with that of the fluctuations. 

Our approach is to use the functional Schrodinger formulation, wherein we specify 
the initial wavefunctional \l/[$(7); t] (or more generally a density matrix p[$(7), $(7); i]), 
and then use the Schrodinger equation to evolve this state in time. We can then use 
this state to compute all of the expectation values required in the construction of the 
effective equations of motion for the order parameter of the theory, as well as that 
for the fluctuations. 

One advantage of this approach is that it is truly a dynamical one; we set up 
initial conditions at some time t Q by specifying the initial state and then we follow 



the evolution of the order parameter (j){f) = and of the fluctuations as this state 

evolves in time. Another advantage is that it allows for departures from equilibrium. 
Thus, issues concerning the restoration of symmetries in the early universe can be 
addressed in a much more general setting. 

Our analysis applies quite generally to any arbitrary spatially flat FRW cosmology. 
We also determine the time evolution of the initial (Boltzmann) distribution functions, 
relate the time evolution to "squeezed states" and perform a numerical integration in 
the case of free fields for de Sitter and radiation dominated cosmologies. We expect 
to provide a quantitative analysis of the evolution of the order parameter and the 
dynamics of phase transitions for interacting fields in a forthcoming article ||. 

The initial state we pick for the field t) is that corresponding to a thermal 
density matrix centered at 4>(t). It is then useful to try to understand the high 
temperature limit of our calculations. We are able to compute both the leading and 
subleading terms in the high T expansion of (0 2 (t)). From this we show that the high 
T expansion cannot be valid for all time, but breaks down in the large time limit. We 
then compute the limit of long times (t — > oo). during the phase transition. 

The time evolution of the Boltzmann distribution functions (initially the thermal 
equilibrium distribution functions) is obtained in ref.0. It is pointed out that to one- 
loop order and also in the Hartree approximation, the time evolved density matrix 
describes quantum "squeezed" states and time evolution corresponds to a Bogoliubov 
transformation. 

To illustrate the departure of equilibrium, we have studied numerically in ref.]7| 
the case of a free massive scalar field in de Sitter and radiation dominated cosmologies. 
It was found that a suitably defined coarse-grained non-equilibrium entropy (per k 
mode) grows linearly with time in the de Sitter case but it is not a monotonically 
increasing function of time in the radiation dominated case. This result may cast 
some doubt on the applicability of this definition of the non-equilibrium entropy. 
There still remain some (open) fundamental questions regarding the connection of 
this entropy and the thermodynamic entropy of the universe, in particular whether 
the amount of entropy produced is consistent with the current bounds. 

This work sets the stage for a numerical study of the dynamics of phase transitions 
in cosmology fully incorporating the non-equilibrium aspects in the evolution of the 
order parameter and which at the same time can account for the dynamics of the 
fluctuations which will necesarily become very important during the phase transition. 

We have applied similar methods to investigate the formation of disordered chiral 
condensates in high energy collisions |J. 

2 Evolution Equations 

We consider the inflaton scalar field <&(x,t) in the spatially flat FRW cosmology (0) 
with action given by 




t)f 
2 a(ty 



V(®(x,t)) 



(3) 




(4) 



with 1Z the Ricci scalar and £ a non- minimal coupling constant. 

In the Schrodinger representation (at an arbitrary fixed time t), the Hamiltonian 
becomes 

Since we consider a 'thermal ensemble', we work with a (functional) density matrix 
p with matrix elements in the Schrodinger representation p[$(7), $(T); t]. We will 
assume that the density matrix obeys the functional Liouville equation 

zhf t =[H(t),p] (7) 

Normalizing the density matrix such that Trp = 1, we define as 'order parameter' 

<f>(t) = ^ / rf 3 x($(f , t)) = ^J d 3 x Tr [p(f)$(£ )] (8) 

where f2 is the co moving volume. 

The evolution equations for the order parameter are as follows: 

d 2 (j){t) ^ 3 tZa(t) , 1 f ,5V ($) 



rft 2 a(t) tZt fi7 N 5$(x) 
It is now convenient to write the field in the Schrodinger picture as 

= (j)(t)+ri(x,t) (10) 
M*,*)> = (11) 

Expanding the right hand side of @ we find the effective equation of motion for 
the order parameter: 

^ + ^m +nm+ yxj^ m+ .... a (12) 

where primes stand for derivatives with respect to <fi. 

This equation of motion is clearly very different from the one obtained by using 
the effective potential. It may be easily seen (by writing the effective action as the 
classical action plus the logarithm of the determinant of the quadratic fluctuation 
operator) that this is the equation of motion obtained by the variation of the one- 
loop effective action. 

The static effective potential is clearly not the appropriate quantity to use to 
describe scalar field dynamics in an expanding universe. Although there may be 
some time regime in which the time evolution is slow and fluctuations rather small, 
this will certainly not be the case at the onset of a phase transition. As the phase 
transition takes place, fluctuations become dominant and grow in time signaling the 
onset of long range correlations^, [|. 



3 Hartree equations 



Motivated by our previous studies in Minkowski space |5|, U which showed that the 
growth of correlation and enhancement of fluctuations during a phase transition may 



not be described perturbatively, we proceeded to obtaining the equations of motion in 
a Hartree approximation J7J. This approximation is non-perturbative in the sense that 
it sums up infinitely many diagrams of the cactus or bubble-type|3|]. The Hartree ap- 
proximation becomes exact in the N — > oo limit of an O(N) vector theory. It at least 
provides a consistent non-perturbative framework in which correlations and fluctua- 
tions can be studied. Corrections to the Hartree approximation can be systematically 
implemented in a consistent fashion. 

The Hartree self-consistent approximation is implemented as follows. We decom- 
pose the field as in fllCf) , using a potential as in (p. We find that the Hamiltonian 
becomes 



+ ^ A ^^ 3 + i! A ^)} ( 13 ) 

The Hartree approximation is obtained by assuming the factorization 

V 3 {x,t) -> 3(r ? 2 (x,t))r/(f,t) (14) 

V \x,t) -> 6( V 2 (x,t))r ] 2 (x,t)-3(r ] 2 (x,t)) 2 (15) 



where (■ • •) is the average using the time evolved density matrix. This average will 
be determined self-consistently (see below). Translational invariance shows that 
(r] 2 (x,t)) can only be a function of time. This approximation makes the Hamilto- 
nian quadratic at the expense of a self-consistent condition. In the time independent 
(Minkowski) case this approximation sums up all the daisy (or cactus) diagrams and 
leads to the self-consistent gap equation 0. In this approximation the Hamiltonian 
becomes 

H = Q a 3 (t) V(0) + (16) 

i A {-Sit) w + f + a3(< » ( vw w " + 5 v,2) <^" 2 ) } 

v(<» = v(0) - Imti 2 } 2 , v«(«) = v'W) + \\<l> iv 2 ) , v< 2, (,» = y"(<» + ±\{ r, 2 ) . 

It is convenient to introduce the discrete Fourier transform of the fields in the 
comoving frame as 

^(£,t)= * £ % (*) e -^ (17) 

The Hamiltonian ( |T6"D then takes the form 
H = Q a 3 (t) V((f)(t)) + 

where the time dependent frequencies (u}% (t) ) and the linear term in rj have the values 

<4(t) = a{t)k 2 + a 3 {t)V i2) {<P(t)) (19) 
V?\<j>(t)) = V«(0(t))V^% o (20) 



We propose the following Gaussian ansatz for the functional density matrix elements 
in the Schrodinger representation 



Vk(t) 



JjA4(t) exp 

k 



Afl \(t) V - k (t) + ^Ut)v- k (t) + ^Vk(t)V-k(t) 



2h 



2h 



+t Kk{t) {r)-k{t) - fj-k(t)) 
h 



(21) 



where <fi(t) = (&{x)) and 7Tfe(i) is the Fourier transform of (U(x)). This form of the 
density matrix is dictated by the hermiticity condition p^[$,$,t] = p*[$,$,t]; as a 
result of this, B k (t) is real. The kernel Bk(t) determines the amount of "mixing" in 
the density matrix, since if B k = 0, the density matrix corresponds to a pure state 
because it is a wave functional times its complex conjugate. 

In order to solve for the time evolution of the density matrix (0) we need to specify 
the density matrix at some initial time t . It is at this point that we have to assume 
some physically motivated initial condition. We believe that this is a subtle point 
that has not received proper consideration in the literature. A system in thermal 
equilibrium has time-independent ensemble averages (as the evolution Hamiltonian 
commutes with the density matrix) and there is no memory of any initial state. 
However, in a time dependent background, the density matrix will evolve in time, 
departing from the equilibrium state and correlation functions or expectation values 
may depend on details of the initial state. 

We will assume that at early times the initial density matrix is thermal for the 
modes that diagonalize the Hamiltonian at t D (we call these the adiabatic modes). 
The effective temperature for these modes is ksTo = l/(3 . It is only in this initial 
state that the notion of "temperature" is meaningful. As the system departs from 
equilibrium one cannot define a thermodynamic temperature. Thus in this case the 
"temperature" refers to the temperature defined in the initial state. 

Inserting the Gaussian Ansatz into the Liouville equation yields upon using 
the hamiltonian (|16|) the following equations for the functions Ak(t), Bk(t) and J\fk(t): 



M 
iA k 

iB h 



2aW) {Ak ~ Al) 



a (ty 



u 2 k(t) 



-^——{A k — A* k ) 
a 3 {t) 



(22) 
(23) 
(24) 



The equation for B k {t) reflects the fact that a pure state B k = remains pure under 
time evolution. 

Writing A k in terms of its real and imaginary components A k (t) = A^ k {t) +iAi k {t) 
(and because B k is real) we find that 



B k {t) B k {t ) 



A Rk {t) A Rk {t ) 
and that the time evolution is unitary (as it should be), that is 

N k (t) 



{A Rk {t) + B k {t)) 



constant 



(25) 



(26) 



For the chosen thermal initial conditions at t = tn, we have 



A k (t ) = A* k (t ) = W k (t ) a 3 (t ) coth\j3 hW k (t )} 
m(t )a 3 (t ) 



where 



B k (t Q 

m(t ) 



sinh \j3 KW k (to)] 

m(toW(t ) 



k 2 + m 2 (T ) 



tanh 



PohW k (t ) 



m\T ) 



v {2) m )) 



(27) 



(2* 



(29) 



a{t ) a? (t ) 

The initial density matrix is normalized such that Trp(t D ) = 1. Since time evolution 
is unitary such a normalization will be constant in time. For T Q = the density 
matrix describes a pure state since B k = 0. 

It is convenient to introduce the complex function: 



A k {t) = tanh %hW k {t )} ReA k (t) + i lmA k (t) 

Then A k {t) obeys a Riccati equation 

1 



iA k {t) 



a 3 (t) 



A h \t) -ui{t)a\t) 



with the initial conditions: 



A k (t ) = W k {t )a 3 {t ) 
Eq.(|3TD can be easily linearized by introducing the functions <f k (t) as 

Mt) = _ ia 3 (i) || 

The equal time two-point function thus becomes 

h f d 3 k 



(30) 

(31) 
(32) 

(33) 



<p£(t) | 2 coth[/? ^>V fc (t )/2] 



(34) 



which leads to the following set of self-consistent time dependent Hartree equations: 



<f) + 3-<f) + V'{(f)) + X(f) n 
a 2 



h f d 3 k I (p%(t) 



(27T) 



coth[(3 hm{t )/2] =0 (35) 



<Pk(to. 



h f d*k 1 yf (t) 
(2vr) 3 2 

1 



Ja?(t )W k (t 



coth[p o hW k (t )/2] 
, <fik(t) \ t = i 



tf® = 



m(to) 



(36) 
(37) 



That is, we are faced with two complicated coupled non-linear and non-local equations 
for (f)(t) and (fk(t). 



The Hartree self-consistent equations ( pop need regular izat ion and renormalization 
since the momentum integrals of the mode functions (^) diverge. 

Because the Bose-Einstein distribution functions are exponentially suppressed at 
large momenta, the finite temperature contribution will always be convergent and we 
need only address the zero temperature contribution. 

A WKB analysis of eq. (|35|) provides the large k behaviour of the mode functions 

<Pk(t) E- 



2a 2 (t)k 4P 



+0(l/k 4 )+- 



Introducing an upper momentum cut-off A we obtain 



(3? 



(V 2 (x,t)) 



h A 2 h 



lTT 2 a 2 (t) 8vr 



-4 



+ finite 
(39) 

where we have introduced a renormalization point K, and the finite part depends on 
time, temperature and K. 

There are several physically important features of the divergent structure obtained 
above. First, the quadratically divergent term reflects the fact that the physical 
momentum cut-off is being red-shifted by the expansion. (This term will not appear 
in dimensional regularization). 

Secondly, the logarithmic divergence contains a term that reflects the initial condi- 
tion (the derivative of the expansion factor at the initial time t a ) . The initial condition 
breaks any remnant symmetry. For example, in de Sitter space there is still invari- 
ance under the de Sitter group, but this is also broken by the initial condition at an 
arbitrary time t D .Thus this term is not forbidden, and its appearance does not come 
as a surprise. As a consequence of this term, we need a time dependent term in the 
bare mass proportional to l/a 2 (t). 

We are now in a position to present the renormalization prescription within the 
Hartree approximation. The differential equation for the mode functions ( j36[) must 
be finite. Thus the renormalization conditions are obtained from 

Afl j2/i \ , * ^ , Xr 



m, 



2 w ^ + W+^ + ir ( ™ (40) 

where the subscripts B, R refer to bare and renormalized quantities respectively and 
(i] 2 ) b is read from (|39|) 



(V 2 )s 

h 



h A 2 



S7T- 



a 2 (t) 



(41) 



In — 



K 



a 2 (t ) 
a 2 {t) 



n 



A 



+ finite 



Using the renormalization conditions ( flop we obtain 



m 2 B (t) + 



X B h A 2 
16tt 2 a?{t) 
Xr 



+ ^4 In 



16tt 2 



a 2 (t ) 
a 2 {t) 



m 



R 



1 + ^4 In 

16tt 2 



Zb = Zr + 



X B h 



(v 2 ) 



R 



167T 2 

Ir + J 



ln(A 

K 



Zr-\ 



(42) 
(43) 

(44) 
(45) 



where 



J 



d 3 k 

6{k - K) 

+ r d 3 k 
n 



*f(*) 



K 

~6 



a 2 (t ] 
a 2 (t) 

vf (t) I 2 



2ka 2 (t) 



A 



^ii / 2\ 



(2tt)3 exp P hW k (t ) 



(46) 
(47) 



The conformal coupling £ = 1/6 is a /ixec? point under renormalization |10[| . In 
dimensional regularization the terms involving A 2 are absent and In A is replaced by a 
simple pole at the physical dimension. Even in such a regularization scheme, however, 
a time dependent bare mass is needed. The presence of this new renormalization 
allows us to introduce a new renormalized mass term of the form 



a 2 (t) 

This counterterm may be interpreted as a squared mass red-shifted by the expansion 
of the universe. However, we shall set £ = for simplicity. 

Notice that there is a weak cut-off dependence on the effective equation of motion 
for the order parameter. 

For fixed Xr, as the cutoff A — > oo 



(47T) 



Mr 



IB 



+ 



In A 



m 



lit) 



o?(t) 



A 5 



+ a 2 (t 



+ 



In A 



(4$ 



This approach to 0~ of the bare coupling as the cutoff is removed translates into an 
instability in the bare theory. This is a consequence of the fact that the N-component 
$ 4 theory for — > oo is asymptotically free (see ref.JTl]]), which is not relieved in 
curved space-time. Clearly this theory is sensible only as a low-energy cut-off effective 
theory, and it is in this restricted sense that we will ignore the weak cut-off dependence 
and neglect the term proportional to the bare coupling in (|49D . 

The renormalized self-consistent Hartree equations thus become after letting A = 

oo: 



+ 3 - a j> + m\<\> + i R n <p + y 
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(49) 
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(v 2 ) 
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tf{t) = 



where [r] } r is given by equations (|4 



4 High Temperature Limit 

One of the payoffs of understanding the large-fc behavior of the mode functions (as ob- 
tained in the previous section via the WKB method) is that it permits the evaluation 
of the high temperature limit. We shall perform our analysis of the high temperature 
expansion for the Hartree approximation. 

The finite temperature contribution is determined by the integral 



For large temperature, only momenta k > T contribute. Thus the leading contribu- 
tion is determined by the first term in eq. (|38|) ). We find 



J 



1 



I2h 



k B T a(t ) 
a(t) 



[1 + 0(1/T ) 



(51) 



Thus we see that the leading high temperature behavior results in an effective 
time dependent temperature 



Teff(t) = T Q 



a{t ) 
a(t) 



This expression corresponds to what would be obtained for an adiabatic (isen- 
tropic) expansion for radiation (massless particles) evolving in the cosmological back- 
ground. 

This behavior only appears at leading order in the high temperature expansion. 
There are subleading terms that we computed subsequently J/J. To avoid cluttering 
of notation, we will set kg = h = 1 in what follows. 

We define 



(52) 



and we will assume that m 2 {T ) -C T 2 . Since we are interested in the description of 
a phase transition, we will write 
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Agg . 
24 : 



T 2 > 



Thus, to leading order in T Q 



m 2 {T ) 



— (T 2 — T 2 ) 

V o c) 



(53) 



(54) 



Our high temperature expansion will assume fixed m{T ) and m(T )/T 1. 
It becomes convenient to define the variable 



A- 2 



x 



m 2 (T 
T 2 a 2 ( to ) ' T 2 
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<pg(t) | 2 = ^ H (a(t )^x 2 T 2 -m 2 (T );t) 



Recall from our WKB analysis that the leading behavior for k — > oo is (see equation 
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2 2a 2 (t)k 

adding and subtracting this leading term in the integral J and performing the above 
change of variables, we have 

J = J1 + J2 
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(55) 



We now must study the high temperature expansion of J\. We restrict ourselves to 
the determination of the linear and logarithmic dependence on T D . For this purpose, 
it becomes convenient to introduce vet another change of variables x = and 

-L o 

use the fact that in the limit T Q 3> m(T ), 
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e T ° ' 
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This yields the following linear and logarithmic terms in T Q : 
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(56) 



Note that the above integral is finite. 

The logarithmic contribution is obtained by keeping the 0(1/ k 3 ) in the large 
momentum expansion of | fk{t) \ 2 given by equation (|38|) (in terms of the new 
variable z). We obtain after some straightforward algebra: 
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(57) 



That is, in the limit T Q » m(T ), J x = J llin + J llog + 0((T ) ). 

Comparing the 0(T 2 , T Q , lnT ) contributions it becomes clear that they have very 
different time dependences through the scale factor a(t). Thus the high tempera- 
ture expansion as presented will not remain accurate at large times since the term 
quadratic in T Q may become of the same order or smaller than the linear or loga- 
rithmic terms. The high temperature expansion and the long time limit are thus not 
interchangeable, and any high temperature expansion is thus bound to be valid only 
within some time regime that depends on the initial value of the temperature and the 
initial conditions. 

As an illustration, we calculated J mn explicitly in the case of de Sitter space [171. 



5 Large Time Limit 



The expansion factor a(t) tends to infinity in the limit t — > oo. As we have just seen, 
the high temperature expansion breaks down in this limit. Physical momenta are 
given by 

/ k 



a(t) 



(hi 



For large a(t), only comoving momenta k = I a(t) — > oo will be relevant. Thus, we can 
again use the WKB method to evaluate the mode functions ^^(^ i n this regime. Let 



us consider for example the de Sitter universe (a(t) = a e ). We find from eq. 
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for m/if << 1. Using the asymptotic behaviour ([59]) in eqs.(|45| ) - (|47|) leads to 
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where 
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(62) 
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It should be noticed that the integral in Y(T) is the mean value of | for a free 



Bose gas at temperature T . Notice that both YiT) and Z( m(Tn) ) depend on the 
renormalization point K and that Z{ 



K 



m(T ) ■ 



m(T ) ■ 

is positive for all K . 



Assuming that (f>(t) — > 0oo (a constant) for £ — ■> oo, the Hartree equations ([49] 
yield for t — > oo, 



> 2 

oo 



A , 



m 2 + 12iP£ r 



-H 2 Z{ 



m(T Q ) 



(63) 



Since Z( m ^ T ^ ) > and we consider X r « 1, this equation has a solution provided 



m 2 + 12# 2 £ f . < 



(64) 



The physical interpretation is as follows, in the symmetry breaking case where eq. ( |5l| ) 
holds, the order parameter <f>(t) for large times t (when the universe cools down fast) 
tends to the non-zero value (j)^ given by eq.(^) independent of the initial value <fi(to). 
The equal time two-point field correlator is given in the Hartree approximation 

by 



A(x,t) = (r)(x,t)r){0,t)) 



d 3 k 



Ak.x I , .H 



<pZ(t) | 2 coth[/U>V fc (t )/2], (65) 



For long times [ a(t) » 1 ] we find using the previous WKB results for | ^{t) \ 2 



A(x, t) 



<*(*)-«> (H\ 2 A(r) 



+ 



B(r) 



where 



(66) 



2tt/ r [2Trra(t)} 2 

71 f°° 

A(r) = / dr Ko(mr) , B(r) = mrKi{mr) , r = Ixl (67) 

2m J r 

We find in terms of the physical (redshifted) length R = ra(t) when t —>■ oo and R is 
kept fixed: 



A{x,t) 



a(i)— >oc 1 



+ 



H 



(2nR) 2 ' "2tt 

Notice that m = m(To) fixes the scale of the correlations. 



[Ht - \og{2mR) + 1 - 7] + 0{\oga{t)/[a{t)} 2 ) (68) 
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